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ABSTRACT 



Optically selected clusters of galaxies display a relation between their optical mass estimates and their X-ray luminosities 
Lx that has a large scatter. A substantial fraction of optically selected clusters have Lx estimates or upper limits 
significantly below the values expected from the Lx-mass relation established for X-ray selected clusters, i.e., these 
clusters are X-ray underluminous for their mass. We attempt to confirm or falsify the X-ray underluminous nature of 
two clusters, Abell 315 and Abell 1456, by using weak gravitational lensing as a third and independent measure of the 
clusters' masses. We obtained optical wide-field imaging data and selected background galaxies using their colors and 
measured the shear exerted by the tidal field of the foreground galaxy clusters. We then fitted parametrized models 
to our shear catalogs. After accounting for projections of large-scale structure and halo triaxiality, we find that A 315 
is significantly X-ray underluminous for its mass, while no significant lensing signal was detected for A 1456. We re- 
evaluate earlier kinematic and X-ray analyses of these two clusters and discuss the nature of the X-ray underluminous 
cluster A 315 and why A 1456 was probably erroneously identified as being X-ray underluminous. 

Key words, gravitational lensing - galaxies: clusters: individual: A 315 - galaxies: clusters: individual: A 1456 



1. Introduction 

Clusters o f galaxies c an be detected by vari ous methods: op- 
tical (e.g., [Ga3|2006j), X -ray emission (e.g., Bohrin ger et al] 
20011). weak g ravitational lensing (e.g., ISchneiderT 11996: 
Dietrich et al] 120071). and the Sunyaev-Zeldovich effect 



McNamara et al.ll200lt iDonahue et alJl2002l iGilbank et al 
2004 iLubin et all 12004: iBarkhouse et al J 120061 : iFang et al 



2007tlPopesso et al.ll2007al P07 hereinafter). These clusters 



(e.g., Staniszcwski ct al. 2008). Differences in the selection 
method could potentially lead to biases when determin- 
ing the cluster mass function. Optical selection is generally 
more affected by projection effects than X-ray selection, 
although projection effects can be minimized b y selecting 
clust er galaxies on the basis of their colors (e.g., Em et all 
20021) . On the other hand, X-ray selection requires that the 
intra-cluster gas has been heated to a detectable level, and 
theoretical predictions show that there is a non-negligible 
fraction of unvir ialized, relatively massive clusters with no 
X-ray emission (<Weinberg fc Kamionkowskj||2002D . While 
such systems are mostly expected at high redshifts, they 
may be present at all epochs. 

Most cluster mass functi ons have been obtained from 
X-ray selected samples (e.g.. iReiprich fc Bohringeil 12002). 
X-ray selection is generally considered to be well under- 
stood, almost pure and complete in mass. However, sev- 
eral investigations have questioned the completeness of 
X-ray selected cluster samples. These investigations have 
shown that there is a population of optically selected 
clusters that deviate from the relation between the X- 
ray luminosity, Lx, and virial mass (or an optical mass 
proxy such as richness and velocity dis persion) established 
for X-ray selected cluster samples (jBower et al.1 Il997t 



are underluminous in X-ray for their masses. P07, in partic- 
ular, identified X-ray underluminous clusters in the RASS- 
SDSS survey among Abell clusters, which were therefore 
named "Abell X-ray Underluminous" (AXU). The virial 
masses of these AXU clusters were determined from the red- 
shifts of their member galaxies, using data from the Sloan 
Digital Sky Survey. X-ray luminosities or upper limits to 
X-ray luminosities were derived from the ROSAT All-Sky 
Survey. 

What is the nature of these AXU cluster s? By analyzing 
the velocity distribution of their members, lP07l suggested 
they are systems in the process of formation, . However, 
it is possible that at least part of the scatter in the Lx~ 
mass relation is not intrinsic, but originates in erroneous 
estimates of either the cluster Lx or its mass. To inves- 
tigate this point, we obtained optical wide-field observa- 
tions of two AXU clusters with the aim of measuring these 
clusters' masses by means of weak lensing, and deriving a 
third, independent mass estimate. We describe these data 
and their reduction in Sect. [5] Section [3] contains the weak 
lensing analysis, which we compare with revised kinematic 
and X-ray estimates in Sects. 0] and We revisit the optical 
cluster luminosity and richness in Sect. [6] and present our 
conclusions in Sect. [71 

Throughout this paper, we use a standard ACDM 
cosmology with fl m = 0.3, Q\ — 0.7, and Hq = 



2 
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70 h km s 1 Mpc 1 . Confidence intervals correspond to the 
68% confidence level. 



2. Data 



Table 1. Coefficients of the lensfit prior function. 



Abell 315, a galaxy cluster at z = 0.174 (|P07t) . was observed 
with the Wide-Field Imager (WFI) at the ESO/MPG-2.2m 
telescope in B-, V-, and R-band. WFI is a focal-reducer type 
camera with a 4 x2 mosaic of 2k x 4k CCDs with a filling 
factor of 95.9% (|Baade et al.lll999f ). Its field-of-view (FOV) 
is 34' x 33', resulting in a pixel scale of 0'.'238. At the cluster 
redshift, 1' corresponds to a physical scale of 177 h^ 1 Mpc. 
The observations were carried out in service mode during 
the nights from Nov. 5 to 12, 2007 in dark and clear sky 
conditions. The total exposure times for the three bands 
were 2880 s, 5890 s, and 5500 s, respectively. 

Abell 1456, which is at a redshift of z = 0.135 (|P07f ). 
was also observed with WFI in the same passbands. For 
A 1456, 1' corresponds to 144 /i -1 Mpc. These observations 
were performed during the nights of May 9 and 10, 2008 in 
dark and clear conditions. The total exposure times for the 
A 1456 field were 1500 s, 5750 s, and 5500 s for the B-, V-, 
and R-band, respectively. 

All data were process ed with the GaBoDS/THELI 
pipeline (jErben et al.ll2005h . Because no photometric stan- 
dard stars were observed in the nights that our data were 
taken, the transformation equations of iKoch et al.l ([2004) 
were used to calibrate the V- and R-band data using SDSS 
magnitudes of objects in the same fields. The zero points 
of the B-band dat a were fixed b y matching the expected 
stellar colors of the iPicklesI (|1998l ) stellar library in a color- 
color diagram to the observed colors of stars in the fields. 
The color s of galaxies were corr ected for Galactic extinction 
using the Schlcgcl et al.l <| 199 8?) extinction maps. The effec- 



tive seeing of the coadded images R-band images is 0'.'74 
(A 315) and l'.'O (A 1456). 

The observations were divided into observing blocks 
(OBs) of 5 dithered exposures. Additional offsets between 
OBs ensured that the sky coverage of the observations was 
more homogeneous. Every region in the final V- and R-band 
images was covered by at least six exposures. In the shorter 
B-band observations, which used only one OB, every lo- 
cation was covered by at least three exposures. The WFI 
point-spread-function (PSF) is smooth across chip gaps and 
has only slow spatial variations, so that the coaddition of 
dithered exposures does not lead to discontinues in the PSF 
within the final image. 

2.1. L ensing ca ta logs 

The Bayesian shape fittin g method /ensfid] (jMiller et alj 
[20071 IKitching et all [2008) was used to estimate the shear 
signal imprinted on the shapes of background galaxies. 
Stars were pre-selected in a magnitude-flux-radius dia- 
gram with the additional requirement that SExtractor's 
CLASS_STAR > 0.95. These stars were modeled with 
the elliptical Gauss-Laguerre met hod ([Bernstein fc Jarvisl 
2002; iNakajima &: Bernstein! l2~007h and a two dimensional 
4th order polynomial was fitted to the Gauss-Laguerre co- 
efficients to obtain a PSF model over the entire WFI field- 
of-view. Lensfit - in the implementation used here - sub- 
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STEP1 


A 315 


A 1456 
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0.029 


0.042 


0.068 
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0.45 


0.44 


0.44 
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42.7 


32.6 


30.5 



divides the image area into regular grid cells, inside which 
the PSF is assumed to be constant. Here the side length of 
the cells was chosen to be 100 pixels or 0f4. For each cell, 
an image of the PSF was reconstructed from the shapelet 
coefficients at the cell's center. These PSF images were then 
used as a PSF description for lensfit. 

Th e Bayesian prior function proposed by IKitching et all 

(gmi, 



V(ex, e%) — A cos 



exp 



2|e| 



B(l 



(1) 



1 If the reader wishes to use the lensfit code, its authors can 
be contacted at http://www.physics.ox.ac.uk/lensfit/ 



where B, C, and D are free parameters and \e\ — \J e\ + ef, 
was used in computing the posterior ellipticity probabil- 
ity distribution. Because of the many degeneracies of this 
function, a brute-force minimization was the only viable op- 
tion in determining its coefficients. The co efficients found 
here a re very similar to those ob tained by | Kitching et alj 
([2008? ) for the STEP1 simulations ([Hevmans et al.ll2006l . see 
Table [1] for a comparison). 

Shear estimat es for individual gala xies were obtained 
using Eq. (20) of IKitching et~aTl (j2008h . When measuring 
the shear we must - as in Eq. ([T]) - assume a prior that 
contains zero shear because we cannot predict the variation 
in shear across the FOV. Since after lensing the ellipticity 
distribution of galaxies is not centered on zero anymore, we 
need to apply an additional weighting factor to correct for 
the assumption of zero shear. The shear sensitivity dei/dgi 
describes how the measured ellipticity e of a galaxy depends 
on the reduced shear g. Very low signal-to-noise ratio (SNR) 
objects have sensitivity values close to zero. Together with 
high ellipticity measurements, this can lead to unphysically 
high shear values. Therefore, shear estimates with shear 
sensitivity of dei/dgi < 0.3 or ellipticity |e| > 1.1 were 
rejected. We note that both types of rejection to low SNR 
objects affected fewer than 50 galaxies, and hence did not 
bias our lensing estimates. 

Because this is the first application of lensfit to real 
data, we compared the shear estima tes obtained i n this 
way with those comp uted using t he Erben et~aT1 {2001) 
implementation of the iKaiser et al.l <| 19951 KSB) shape es- 
timation algorithm. The STEP1 bias parameters q (non- 
linear response), m (calibration bias), and c (constant off- 
set) were computed where the lensfit shear estimates were 
set to 7* rue . We found that q and c are consistent with 
zero and that m = 0.152 ± 0.029. This is consistent with 
the bia s parameter values of the MH implementation of 
KSB in lHevmans et alj (|2006f ). which is, save a few details, 
identical to the KSB implementation used here. This result 
suggests that, apart from errors in the PSF model common 
to both methods, the lensfit shear estimates are very close 
to the true shear values. 

Background galaxies were selected based on their mag- 
nitude and colors. Galaxies with R > 23 mag were consid- 
ered to be background galaxies and were included in the 
catalog irrespective of their colors, while brighter galaxies 
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Fig. 1. Left: Mass reconstruction of the A 315 field. The surface mass density contours start at at 0.0464 (corresponding 
to 1.6 x 10 14 hM Q Mpc -2 or 2a above the mean surface mass density at the edge of the field), increasing in steps of la. 
Dashed contours are at the same negative levels. Right: B-mode reconstruction of the A 315 field. Contours are at the 
same levels as in the left panel. 
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Fig. 2. Same as Fig. Q] for the A 1456 field. The la level is at 0.032 (corresponding to 1.25 x 10 14 h Mq Mpc -2 ). The 
color bar on the right runs from —la to 5a in both figures. 



were rejected if their colors matched the colors of cluster 
galaxies, 0.53 < (B — V) < 1.38 and -0.25 < 1.7 x (V - 
R) — (B — V) < 0.3. Finally, only galaxies manually prese- 
lected in a magnitude-flux-radius diagram were kept if their 
flux radius was < 1'.'3. The effectiveness of this color selec- 
tion was tested by plotting the number density of galaxies 
in radial bins from the cluster centers. A small excess of 
galaxies was found only for A 315 within a clustercentric 
distance of D2; beyond this radius, the number density re- 
mained constant, which is indicative of no contamination 
of the background catalogs by cluster galaxies. 



After these selections, 9505 galaxies with shear es- 
timates remained in the lensing catalog of A 315, and 
8030 galaxies were left in the A 1456 catalog, correspond- 
ing to respective number densities of 8.6arcmin~ 2 and 
7.5arcmin -2 , where areas masked due to the occurrence 
of bright stars, reflection rings, and diffraction spikes were 
excluded. 
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3. Weak lensing mass estimates 

Based on the lensing catalogs described in the previous sec- 
tion, surface mass density m aps were computed using t he fi- 
nite field inversion method of lSeitz fc Schneiderl (|2001[ ) with 
a smoothing scale of 2'. These maps are presented in the 
left panels of Figs. Q] and [2] The right panels of these fig- 
ures are B-mode maps obtained from shear catalogs with 
galaxies rotated by 45° to cancel any true shear signal, and 
is a test of systematic residuals in the PSF correction. The 
flatness of the B-mode maps suggests that the PSF models 
used in creating the lensing catalogs are sufficiently accu- 
rate for obtaining cluster mass estimates. The noise levels 
of the reconstructions were estimated from the variance of 
100 maps with randomized galaxy orientations. 

At a confidence level just under 5a, A 315 is clearly 
detected a little SE of the image center. No lensing signal 
is visible at the position of A 1456 at the center of the 
image. The strong peak at the right of the image is probably 
caused by the reflection ring of a bright star on which it is 
centered. However, this peak is also only 45" away from t he 
cluster candidate NSC J120257+040951 (|Gal et al.ll20"03h . 

The weak lensing peak of A 315 is centered on a galaxy 
concentration whose brigh t est ga laxy is 2'.2 from the clus- 
ter position of lAbell et all (|l989h and 4f3 from the cluster 
center adopted in the X-ra y and kinematic analysis of lP07l 
Both the position of IP07I and the brightest cluster galaxy 
(BCG) identified here (02:10:06.46, -01:01:56.4 6), are com- 
patibl e with the positional uncertainties of lAbell et alJ 
1989) but are mutually incompatible. Because the IP071 
position of A 315 was based on a very weak X-ray sig- 
nal, the BCG position was adopted as the cluster center 
for the analysis presented here, a choice that we further 
justify in Sect.lU where the kinematic analysis of both clus- 
ters is re-examined. For A 1456 the weak lensing analysis 
was done taking the BCG identified in Sect. H (12:03: 48.71 , 
+04:20:43.46), located 5f75 from the X-ray center of lP07l 
to be the center. 

Converting the dimensionless lensing quantities into 
physical mass densities requires knowledge of the source 
redshift distribution or at least of the average source red- 
shift. To determine the latter quantity, galaxies were ran- 
domly drawn from the photometric redshift catalog of 
lllbert et all (|2006h to match the magnitude distribution 
of our lensing catalogs, us ing the color transformations 
of iBlanton fc Roweisl (2007). The average redshifts of the 
shear catalogs are ZA315 = 1-07 and zai456 = 1-03. While 
it is formally possible to estimate errors in these mean red- 
shifts by bootstrap resampling from the CFHTLS catalogs, 
the true redshift error will be dominated by cosmic vari- 
ance, which we do not explicitly calculate here. We note 
that for the two low-redshift clusters under investigation, 
a redshift error as high as 0.2 would result in a mass error 
of less than 5% and even in this extreme case would be an 
insignificant contribution to the total error budget of the 
lensing mass estimates. 

To obtain mass estimates, parametrized models - singu- 
lar isothermal spheres (SIS) and NFW pro files - were fitted 
using the maximum likelihood method of ISchneider et al.l 
(|2000fh The best-fit SIS model for A 315 has a velocity 
dispersion of cr S i S = 747tgkms~ 1 . The best-fit NFW 

model for A 315 has M 200 = 2.96±o.75 x 10 u h' 1 M . 
The minimiza tion was carried ou t with a Downhill Simplex 
method (e.g.. IPress et al.l fl992) in which at every vertex 



th e concentrat i on pa rameter was fixed to the prescription 
of iDolag et alJ (|2004l ). Attempts to fit parametrized mod- 
els to A 1456 returned only marginally significant results. 

' and the NFW 
Mp). All model fits 



The best-fit SIS has cr S i S = 418±J£kms 



113 . 



model has M 200 = 6.3+^ x 10 1 
we re centered on the BC G. 

iHoekstral (|200lL I2003D and lDodelsor] (|2004[ ) studied the 
influence of uncorrect ed large-scale str ucture (LSS) on 
cluster mass estimates. iDodelsonl d200l found that LSS 
projections can increase the error in M200 by as much as 
75% when the concentration parameter c and M200 are fit- 
ted simultaneously. Much of this additional error comes 
from the stretching of the error ell ipsis along the c -axis, 
and it was estimated from Fig. 3 of Do delson (2004]) that 
in the case of fixed c, as used in this work, the additional 
err or is at most 50%. 

iCorless fc Kind (|2007f ) investigated the influence of halo 
triaxiality on weak-lensing mass estimates and found that 
very oblate or prolate halos with axis ratio 1:3 decrease 
or increase weak-lensing mass estimates by as much as 
40%, depending on their orientation with respect to the 
line of sight. These are, however, e xtreme cases, which oc- 
cur in less than 1% of all halos (|Kasun fc Evrardl[2005l) 
and a realistic estimation of the contribution to the total 
error budget must t ake th e halo-shape distribution into ac- 
count. [Shaw|eFliI] (|2006l ). for example, fo und that prolate 
halos, i.e., those halos that according to ICorless fe Kind 
(2007) lead to higher additional errors in weak lensing 
measurements, are muc h more common than oblate halos. 
iKasun fe EvrarcH (|2005l ) presented a fitting function for the 
minor-axis as a function of mass and redshift. For a halo 
at the redshift and mass of A 315 the ex p ected c/a = 0.63. 
From Figs. 4 and 5 of ICorless fc Kind (|2007l ). we found 
that in the extreme cases of the major (minor) axis being 
aligned with the line of sight, the weak-lensing mass for 
this axis ratio will be overestimat ed (underestimated) by 
16% (10%). Because the results of ICorless fc Kind (|2007h 
cannot be convolved with the full halo-shape distribution, 
we adopt this combination of mean axis ratio and extreme 
alignment as additional error from halo triaxiality. 

Adding the systematic errors of LSS projections and 
halo triaxiality leads to total errors that are significantly 
higher than those obtained from the model fits alone. For 
A 315 M200 = (2.96±^±^) x lO 14 /^ 1 Af Q was found, 
where the first error is random and the second error is 
sytematic. The marginal detection of A 1456 turns into an 
upper limit to the mass if projection effects of LSS are taken 
into account, M 200 = (6.3±^?±|;|) x 10 13 hr 1 M Q . 



4. Kinematical mass estimates 

The lensing analysis in the previous section established 
that the position of A 315 is differs sig nifica ntly from the 
one used in the kinematical analysis of IP07I . Because the 
choice of centroid position affects the selection of clus- 
ter members, the kinematic ana lysis of SDSS DR5 data 
()Adelman-McCarthv et al.l l2007f) was repeated using the 
B CG as the cluster cen ter. This analysis follows the steps 
of[B iviano et alJ ([2006), as briefly summarized here: 

1. a peak in the redshift distribution is selected; 

2. interlopers are rejected by analyzing the projected 
phase-space diagram, i.e., velocities versus cluster- 
centric distances; 
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Fig. 3. The clusters A 315 (red diamond) and A 1456 
(blue circle) i n com parison with the Lx _ -^200 relation of 
iRvkoff et al.l (|2008l black dots and shaded area). The er- 
ror bars for the mass of A 315 are the sum of random and 
systematic errors. For A1456, we use the la random plus 
systematic upper limit to the mass estimate, since the clus- 
ter is undetected in the weak lensing analysis. 

3. the virial mass and velocity dispersion are computed 
within a given aperture; 

4. the computed mass and velocity dispersion are trans- 
lated into the mass within r2oo- 

The mass estimates are c orrected for the surface-pressure 
term (|The fc Whitelll986l) . In the following, we describe the 
influence of these steps on the kincmatical mass estimates 
of A 315. 

Step 1 requires the choice of an aperture centered on 
the adopted center positi on to select all ga laxies within this 
aper ture. Weighted gap (jBeers et al.lll990f) and the density 
gap (jAdami et al.lll998l ) estimators with a gap value of 4 
were chosen to select the peak position. We note that our 
mass estimates are insensitive to the choice of the gapper 
method and of the gap value. However, the mass value de- 
pends on the initial aperture for the search of the peak. 
Adopting a radius of 2 h~ x Mpc instead of 1 h~ x Mpc causes 
the mass estimate of A 315 to increase by 31%. This is be- 
cause of the inclusion of 3 galaxies in the data-set that are 
excluded when an aperture of 1 h^ 1 Mpc is used. With the 
exception of these 3 galaxies, the identification of interlop- 
ers (step 2) is quite robust, with little room for a different 
selection. 

The mass estimate in step 3 may be computed within an 
aperture that differs from the one used in step 1 by being 
larger. Here one Abell radius, 2.15 ft. -1 Mpc, was adopted. 
The mass within this radius is used to obtain a first esti- 
mate of r2oo , and then to interpolate from the Abell radius 
to r 2 oo using an NFW profile and obtain M 2 oo (step 4). 
Alternatively, instead of using the Abell radius, the mass 
within the aperture corresponding to r2oo can be computed. 
Since the r2oo estimate will change as a result of the aper- 
ture selection, this is done iteratively. The difference be- 
tween the mass estimates of both methods is negligible. 

Rather than using all cluster galaxies, one can use only 
those on the red sequence. In this case, we found that the 
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Fig. 4. Comparison of kinematic and weak-lensing mass es- 
timates of A 315. The filled (blue) circles are virial mass es- 
timates; the filled (red) squares are mass estimates obtained 
from the velocity dispersion. Shown are the estimates for 
the samples of (1) all galaxies plus the 3 outliers (see text 
for details); (2) all galaxies; (3) only red-sequence galaxies; 
and (4) galaxies within r2oo- The numbers on the x-axis 
indicate the sample sizes. The solid (green) line shows the 
weak-lensing mass estimate with error bars (dashed lines) 
for comparison. 



mass estimate decreases by 25%, if the mass estimated is 
obtained within an aperture of 2.15 h~ x Mpc. 

Two mass estimates are reported at the end of each 
analysis. One comes from the virial analysis and the other 
uses the velocity dis persion as a proxy. It is known from nu- 
merical simulations (|Biviano et al.ll2006T ) that the former is 
generally an overestimate and the latter an underestimate, 
so that the two values should bracket the real one. 

Mass errors were computed with the Jackknife tech- 
nique. They range from 23% to 35% for the virial mass, 
but are much higher for the mass derived from the velocity 
dispersion (from 36% to 90%). The various mass estimates 
of A 315 are displayed in Fig. 2] All mass estimates reported 
here are low er th an the value of (6.61 ±2.71) x 10 14 h~ x M & 
reported bv lP07l The difference is due to the different cen- 
tral position adopted here, the BCG rather than the X-ray 
center. This centroid choice also decreases the systematic 
error estimates because the main peak seems to be more 
reliably centered. 

We consider the range in the estimates obtained by 
changing the parameters of the analysis to be systematic 
error, and the average jackknife errors of these estimates 
as random error. The cluster mass is now estimated to be 
(2.7tl[ 1 7 ± L0 ) x lO 14 /^ 1 M Q , where the first error is the 
random error and the second one is systematic. 

According to a lDressler fc Shectmanl (|1988| ) test, there 
is no evidence of subclustcring, but this test is not very 
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Fig. 5. Same as Fig. 2] but for A 1456. The leftmost point 
is for the analysis performed at the X-ray cluster position. 
All other points were computed with the aperture centered 
on the BCG. 



Fig. 6. Velocity histogram of A 1456. Velocities are given 
with respect to the mean velocity. The vertical red line 
denotes the velocity of the BCG. 



reliable for samples of redshift for fewer than 50 galaxies; 
there is also no evidence of a gradient in the R-v plane. 

The kinematic analysis was also repeated for A 1456, 
where it was found that, compared to the A 315 case, the 
mass estimates obtained for A 1456 depend very little on 
details of the membership selection for galaxies in the clus- 
ter field. Variations in (1) the choice of the cluster center, 

(2) the selection of all vs. red-sequence galaxies only, and 

(3) the choice of the limiting aperture for the analysis, lead 
to variations in the virial mass estimates in the range 4.1- 
5.2 x 10 14 /i _1 Mq, and in the velocity dispersion mass es- 
timates in the range 2.5-3.1 x 10 14 h' 1 Mq. These values 
agree with those listed in IP071 but are significantly above 
the mass limit of 1.75 x 10 14 Mq derived from the weak- 
lensing analysis. A comparison of the weak lensing and the 
different kincmatical mass estimates is shown in Fig. [SJ 

There is only marginal indication of subclustering in 
the A 1456 data. There is no gradient in the radius- 
velocity plot, but a marginally significant result of the 
iDressler fe Shectmanl test for the red-sequence galaxies 
(96% probability of sub-clustering). Moreover, the ve- 
locity histogram for all galaxies within an Abell radius 
(2.15 h" 1 Mpc) is not symmetric, with a skewness of —0.62 
that is marginally significant (90-95%). Finally, the BCG 
is offset from the mean velocity by — 231 km s -1 but the 
difference is within 1.6a of the mean velocity. The velocity 
histogram is shown in Fig. [6l 



5. X-ray data analysis 

Based on the weak-lensing and revised kinematical mass 
estimates, we also reanalyzed the X-ray data of A 315 
and A 1456 in the RASS. The cluster temperature was 



calculated from the cl uster mass at r^on u sing the mass- 
temperature relation in I Zhang et al.l (|2008l Table 3). 

The X-ray luminosity, Lx, depends on the square of 
the electron density, n 2 , and the emission coefficient, A, 
Lx oc J n 2 A(T, Z)dV. The luminosity can be derived from 
the X-ray s urface brightness i f the emission coefficient is 
known fe.g- IZhang et al.ll2005f ). The emission coefficient is 
a function of gas temperature and metallicity. It can be 
derived in XSPEC using a spectral model for the cluster. 
Here we used a combined model of "wabs x raymond " 
(jMorrison fc McCammonl 119831 : [Raymond fc Smith! 119771) . 
The "raymond" model describes the cluster thermal com- 
ponent including parameters of the gas temperature, T, gas 
metallicity (assuming Z ~ 0.4Z©), and cluster redshift. The 
"wabs" model is set to the hydrogen column density from 
the L AB survey ( Hartmann fc Burtonlll997l iKalberla et al.l 



I2005D . 2.5 x 10 20 cm 2 for A 315 and 1.65 x 10 20 cm 2 for 
A 1456. 

Both A 315 (ID: 931706) and A 1456 (ID: 931633) were 
observed with the ROSAT PSPC in the RASS. One can 
subtract the background from the X-ray image to obtain 
the cluster image. The cluster surface brightness is the X- 
ray image divided by the exposure map, which was used 
within the virial radius to determine the X-ray luminosity. 

The cluster temperature of A 315 derived from the weak 
lensing mass is 3.33 keV. We used the BCG position to de- 
termine the X-ray luminosity. An attempt to determine 
an X-ray position for A 315 resulted in a non-significant 
detection at (02:10:07.263, -00:59:49.98). The X-ray lu- 
minosity changes by less than 1% if the cluster center is 
shifted to this location. The derived cluster X-ray lumi- 
nosity is 1.17 ± 0.48 x 10 43 /i" 2 crgss" 1 in the 0.1-2.4keV 
band. The X-ray luminosity of the BCG (02:10:06.458, 
—01:01:56.46) is below the background level. We therefore 
used the background level to estimate an upper limit to the 
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X-ray luminosity of the BCG inside a radius of 1.5'. This 
is L x < 4.58 x 10 42 ft" 2 ergs" 1 in the 0.1-2.4 keV band. 

The cluster temperature derived from the kinematical 
mass (M 200 = 2.7 x lO 14 /^ 1 M Q ) is 3.14keV. The derived 
cluster X-ray luminosity is 1.18± 0.48 x 10 43 h~ 2 ergss -1 in 
the 0.1-2.4keV band. The X-ray luminosity varies by not 
more than 1% for the different cluster masses derived in 
Sects. [3] and |U We note that all X-ray luminosities given 
here for A 315 are derived from very marginal detections. 
The dominant error in these luminosities is probably the 
background subtraction and the true luminosity could be 
substantially lower. 

The cluster temperature of A 1456 derived from the 
weak-lensing mass (M 20 o = 0.63to^ x 10 14 hr 1 M ) is 
1.24keV. We used the X-ray cluster center (12:03:45.7, 
+04:15:00) in our analysis. The la error of the X-ray center 
position is 3', i.e., about 2a away from the BCG. The X-ray 
emission at the BCG position is not significantly above the 
background. 

The derived cluster X-ray luminosity is 2.88 ± 0.51 x 
10 43 /i- 2 ergss- 1 in the 0.1-2.4keV band. The X-ray lu- 
minosity at the BCG position (12:03:48.700, +04:20:44.00) 
within a radius of 1.5' is 0.91 ±0.91 x 10 42 h~ 2 ergs -1 in the 
0. 1-2.4 keV band, which is the sum of the cluster luminosity 
in this region and the BCG luminosity. 

Because the weak-lensing mass estimate is compati- 
ble with no mass being present at the cluster location, 
we also used the upper limit to the weak-lensing mass 
(M200 = 1-75 x 10 14 /i -1 Mq) to estimate an upper limit 
of the X-ray luminosity. This gives a cluster tempera- 
ture of 2.37 keV. The derived cluster X-ray luminosity is 
4.43 + 0.65 x 10 43 ft" 2 ergss- 1 in the 0.1-2.4keV band. It 
is this luminosity that we plot in Fig. [3] 

The cluster temperature corresponding to the upper 
value of the kinematical mass (A/200 = 5.18 x 10 14 h~ x M & ) 
is 4.67 keV. The derived cluster X-ray luminosity is 4.32 ± 
0.66 x 10 43 ft- 2 ergss- 1 in the 0.1-2.4keV band. 



6. Optical cluster luminosity and richness 

After establishing that A 315 indeed seems to be X-ray un- 
derluminous for its mass and that A 1456 is is in agreement 
with the normal LX-M200 relation but has an unusually 
high kinematical mass estimate, we now examine whether 
the optical properties of these clusters are in any way excep- 
tional. Specifically, we determined their cluster luminosities 
and richnesses and compared them to the full RASS-SDSS 
cluster catalog. 

The total optical luminosity of a cluster has to be com- 
puted after the subtraction of the foreground and back- 
ground galaxy contamination. We considered two differ- 
ent approaches to the statistical subtraction of the galaxy 
background. We computed the local background number 
counts in an annulus around the cluster and a global back- 
ground number counts from the mean of the magnitude 
number counts determined in five different SDSS sky re- 
gions, randomly chosen, each with an area of 30 deg 2 . In 
our analysis, we show the results obtained using the optical 
luminosity estimated with the second method. The optical 
luminosity is then c omputed following the prescription of 
iPopesso et al.l (|2004[ ). The reader is referred to that paper 
for a detailed discussion of the comparison between optical 
luminosities calculated with different methods. To avoid se- 
lection effects due to the slightly different redshifts of the 




RASS-SDSS 
Abell Clusters 



M 200 /(10»M 3 ) 

Fig. 7. Location of Abell 315 and Abell 1456 in the L op - 
A/200 relation. The optical luminosity is calculated inside 
r2oo and is corrected for contamination due to projection 
effects. The red filled hexagon refers to Abell 315 and the 
red star is the upper limit for Abell 1456. The empty black 
squares and the filled points are the X-ray and opt i cally se - 
lected clusters, respectively, used in lPopesso et al.l (|2007bl ). 
The solid and dashed line are the best-fit lines obtained in 
IPopesso et all (|2007bD 



clusters, the optical luminosity was calculated in the same 
absolute magnitude range for all the clusters. We used an 
absolute magnitude cut of M r < —20, which allowed us 
to sample the cluste r luminosity function down to M* + 2 
(<Popesso et al.ll2005D . 

The cluster richness N ga \ was calculated by summing 
the background-subtracted cluster number counts used to 
calculate L op . This estimate was corrected for projection 
effects in the same w ay as we did for L OD , a ccording to the 
prescription given in lPopesso et alj (|2007bl ). 

The cluster optical luminosity and richness were cal- 
culated within a physical aperture of ^200- According to 
the weak-lensing results, we used r 2 oo = 1.30 ft -1 Mpc for 
A 315. For Abell 1456 we used the upper limit to the virial 
radius, r2oo = 1.10 /i -1 Mpc, corresponding to the weak- 
lensing limit of 1.75 x 10 14 M Q . 

Figures [3 [HI and [5] show the location of A 315 (red filled 
hexagon) and A 1456 (red star) in the L opt — M200, M/L — 
M200 and jVga i — M20 relations, respectively, obtained by 
IPopesso et atl (|2007bD . 



7. Discussion and conclusions 

We have studied two clusters of galaxies that were previ- 
ously reported to be X-ray underluminous on the basis of 
their their kinematically derived masses. To either confirm 
or falsify the X-ray underluminous nature of these clusters, 
we obtained deep multi-color optical imaging to measure 
the cluster masses using the gravitational lens effect. 
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Fig. 8. Location of Abell 315 and Abell 1456 in the 
hhoo/Lop — M200 relation. The symbols are the same as 
in Fig. E] 



1 000 




100 



Fig. 9. Location of Abell 315 and Abell 1456 in the N ga i 
M 2 oo relation. The symbols are the same as in Fig. [7] 



Our lensing analysis of Abell 315 confirmed that this 
system is a massive cluster, with a lensing mass esti- 
mate of M200 = (2.96±J;75±°:5o) x 10 14 /i _1 Af Q . This is 
substantially lower than the kinematical mass es tima te of 
M 200 = 6.61 ± 2.71) x lO 14 ^" 1 M Q reported by GEM We 
could expl ain this discrepancy by the choice of cluster cen- 
ter in lP07l . Shifting the center in the kinematical analysis to 
that of the BCG in this work leads to a revised kinematical 



'200 - (2.7±S3 ±1.0) x IO^^Mq, 
which is in excellent agreement with the weak-lensing mass 
estimate. 

To reliably establish the X-ray underluminous nature 
of A 315, we also reanalyzed the available X-ray data in 
the RASS. The X-ray luminosity L x = (1.18 ± 0.48) x 
10 43 h~ 2 ergs -1 determi ned in this work is some what higher 
than the value given by iPopesso et al.1 (l2007af) but is still 
signifi cantly below the LX-M200 relation of Rvkoff et alj 
(l200l . 

The lensing analysis of Abell 1456 field did not show 
a massive cluster. Our attempts to fit parametric models 
centered on the BCG position returned only marginally sig- 
nificant results, as would be expected from random fluctu- 
ations due to shape noise and the projections of LSS. We 
therefore interpret the lensing measurement as an upper 
limit to the mass of a cluster at this position and redshift, 
-M200 < 1.8 x 10 14 h^ 1 Mq. This is significantly lower than 
the kinematic mass estimate that we find for A 1456, which 
seems to be very robust. 

A possible explanation of the much lower mass deter- 
mined by weak lensing could be that A 1456 is a bimodal 
cluster that the kinematic analysis is unable to model ac- 
curately, given that indications of subclustering are only 
marginal. If the cluster is indeed bimodal, as suggested by 
the velocity histogram, the position of the BCG in-between 
the two velocity peaks suggests that it has been displaced 
from one of the two peaks by an interaction. Hence this 
would imply that the kinematic mass estimate is incorrect. 

The X-ray analysis of A 1456 shows a much clearer de- 
tection than in the case of A 315. By using the upper mass 
limit derived from the weak lensing analysis to determine 
to which radius the X-ray luminosity is meas ured, we can 
lace A 1456 right on the Lx-M 2 oo relation of iRvkoff et al.l 
2008J), and the X-ray luminosity based on the kinematic 
mass places A 1456 slightly above the LX-M200 relation. 

In summary, our results show that A 315 is indeed sig- 
nificantly below the Lx-^200 relation, confirming that it 
is underluminous compared to X-ray selected galaxy clus- 
ters of the same size. We note that the error in the X- 
ray luminosity is probably dominated by the uncertainty 
in the background model and higher than the one quoted 
in Sect. \5\ We could not confirm that A 1456 is an AXU 
cluster. The absence of a significant lensing signal in com- 
bination with the marginal indication of bimodality in the 
velocity histogram make the explanation that this is an un- 
relaxed cluster merger more likely. 

The revised optical luminosity and richness of both clus- 
ters are unremarkable. The optical properties give no indi- 
cation that A 315 is either X-ray underluminous or dynam- 
ically unrelaxed. This is in agreement with |P07| who found 
that the AXU clusters deviate from the X-ray luminosity- 
mass relation but are on the optical luminosity-mass rela- 
tion. 

How unusual is A 315 then? IRvkoff et alj (|2008f) found 
a scatter in luminosity at fixed mass of <J\ n L x \M — 0-4, 
much higher than the scatter in the binned relation shown 
in Fig. [3J The X-ray luminosity of A 315 given its mass 
estimate would imply that the cluster deviates by 6a from 
the Lx-mass relation, but the uncertainty in the A 315 
mass estimates is rather high. However, the kinematical and 
weak-lensing mass estimates are based on different physi- 
cal principles and agree well with each other. A weighted 
average of the two measurements gives a best estimate of 
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M 20 o = (2.85±g:^? ±g;ls) x 10 14 hr 1 M Q , so that the cluster 
deviates by 2a from the Lx-mass relation along the mass 
axis. In all cases, A 315 appears to be an unusual cluster 
compared to the population of X-ray selected clusters that 
define the Lx-mass relation. 

Weak-lensing and optical-mass estimators are both af- 
fected by projections along the line-of-sight (LOS). An ex- 
tendend structure such as a filament along the LOS or the 
superposition of smaller groups, would increase both mass 
estimates. The X-ray luminosity, which depends on the 
square of the density of the intra-cluster medium, would 
not be enhanced very much in this scenario. While with 
present data be used to exclude such scenarios, nothing 
in the data suggests that this is a configuration that we 
observe. A superposition of groups should produce kine- 
matical substructure. With the few redshifts available from 
the SDSS, this cannot be decisively excluded but we found 
no evidence of substructure in Sect. [U A filament is a low 
density enviro nment and as such is dominated by blue late- 
type galaxies (Braglia ct al. 2007|). We found that 12 of the 
14 galaxies inside r2oo brighter than L* are on the clus- 
ter red-sequence. This shows that A 315 is dominated by 
early-type galaxies as one would expect for a high-density 
environment such as a cluster. 

Studying A 315 and similar X-ray underluminous galaxy 
clusters is important to understanding the nature of the 
large scatter in the cluster X-ray luminosity-mass relation. 
If intrinsic, as supported by our weak-lensing mass mea- 
surement of A 315, this scatter implies that X-ray selected 
cluster samples are incomplete samples in terms of mass. 
An analysis of th e sta cked velocity histograms of all RASS- 
SDSS clusters in lP07l showed a leptok urtic distribution, in - 
dicative of systems still in formation (|Woitak et al.| [2005). 
In such a system still in formation, the intra-cluster gas 
itself may not yet have reached its final temperature. 

This scenario could be tested with further spectroscopic 
and X-ray observations. Data from the SDSS do not provide 
enough spectra of cluster members to show deviations from 
Gaussianity in the velocity histogram of individual galaxies. 
A spectroscopic sample of several hundred cluster galaxies, 
which can easily be obtained with modern spectrographs, 
could be used to confirm the leptokurtic velocity distribu- 
tion for an individual cluster and identify infalling substruc- 
ture. Similarly, deeper X-ray observations than available in 
the RASS would allow us to measure the density distri- 
bution of the intra-cluster gas. An improved temperature 
measurement due to higher quality (background) statistics 
and a more self-consistent X-ray analysis that does not rely 
on other methods to determine the truncations, could sig- 
nificantly reduce the error in the X-ray luminosity and pro- 
vide a more definite answer to the question of how unusual 
Abell 315 really is. 

On the theoretical side, SPH simulations of galaxy clus- 
ters could be used to generate mock X-ray observations, and 
semi-analytic models of galaxy formation could be used to 
populate TV-body halos with galaxies. The combination of 
these two methods could provide further insight into the 
observable properties of galaxy clusters that are not yet 
fully virialized. 
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